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Abstract. 

We have analyzed the interaction of a dissipative two level quantum system with high and 
low frequency excitation. The system is continuously and simultaneously irradiated by these 
two waves. If the frequency of the first signal is close to the level separation the response of the 
system exhibits undamped low frequency oscillations whose amplitude has a clear resonance at 
the Rabi frequency with the width being dependent on the damping rates of the system. The 
method can be useful for low frequency Rabi spectroscopy in various physical systems which 
are described by a two level Hamiltonian, such as nuclei spins in NMR, double well quantum 
dots, superconducting flux and charge qubits, etc. As the examples, the application of the 
method to a nuclear spin and to the readout of a flux qubit are briefly discussed. 



It is well known that under resonant irradiation a quantum two level system (TLS) can 
undergo coherent (Rabi) oscillations. The frequency of these oscillations is proportional to 
the amplitude of the resonant field [1] and is much lower than the gap frequency of TLS. The 
effect is widely used in molecular beam spectroscopy [2] , and in quantum optics [3] . 

During the last several years it has been proven experimentally that Rabi spectroscopy can 
serve as a valuable tool for the determination of relaxation times in solid state quantum me- 
chanical two-level systems, qubits, to be used for quantum information processing [4]. These 
systems normally are strongly coupled to the environment, which results in the fast damping 
of Rabi oscillations. It prevents the use of conventional continuous measurements schemes for 
their detection, though the special schemes for the detection of coherent oscillations through 
a weak continuous measurement of a TLS were proposed in [5-7]. That is why Rabi oscilla- 
tions are measured with the pulse technique through the statistic of switching events of the 
occupation probability between two energy levels with excitation and read out being taken at 
the gap frequency of TLS, which normally, lies in GHz range [8-11]. The main drawback of 
this technique is that it requires sophisticated high frequency readout electronics. 

In this work we propose a new experimental method for a Rabi spectroscopy of a TLS, 
when readout electronics is continuously swept across the low (compare to the gap) Rabi 
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resonance. The method is rather general and can be applied to a great variety of two level 
systems. 

We consider a TLS which is irradiated continuously by two external sources. The first 
with a frequency luq, which is close to the energy gap between the two levels, excites the low 
frequency Rabi oscillations. Normally, Rabi oscillations are damped out with a rate, which is 
dependent on how strongly the system is coupled to the environment. However, if a second low 
frequency source is applied simultaneously to TLS it responds with persistent low frequency 
oscillations. The amplitude of these low frequency oscillations has a resonance at the Rabi 
frequency with the width being dependent on the damping rates of the system. Note, that this 
approach has a well known classical analog. Indeed, a damping rate of a classical oscillator 
can be easily obtained from its amplitude-frequency characteristics. 

We start with a Hamiltonian of a driven TLS , which is subjected to both high and low 
frequency excitation: 
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Here the first two terms describe an isolated TLS, which can model a great variety of situations 
in physics and chemistry: from a spin 1/2 particle in a magnetic field to superconducting flux 
and charge qubits [4], [12]. In order to be exact we consider the first two terms in to 
describe a double-well system where only the ground states of the two wells are occupied, 
with A being the energy splitting of a symmetric (e = 0) TLS due to quantum tunnelling 
between two wells. The quantity e is the bias, the external energy parameter which makes 
the system asymmetric. The last two terms in Q describe the interaction with external 
time-dependent high frequency, F, and low frequency, G fields which modulate the energy 
asymmetry between the two wells. 

Hamiltonian is written in the localized state basis, i.e., in the basis of states localized 
in each well. In terms of the eigenstates basis, which we denote by upper case subscripts for 
the Pauli matrices ax , ay , a z , Hamiltonian Q reads: 
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where A £ = VA 2 + e 2 is the gap between two energy states. 

The inclusion of the dissipative environment in Hamiltonian (JJJ, results in the Bloch- 
Redfield equations for the matrix operators ax, ay, <Jz, [13], [14]. For weak driving (F <C A), 
and for weak coupling of the TLS to the bath these equations can be approximated by Bloch- 
type equations [15]: 



(&z) = (2/coswoi + 2. 9 (£)) (ay) - T z ((a z ) - Z ) , 
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where / = AF/hA E , g(t) = AG(t)/fiA e , and Z = - tanh (A e /k B T) is the equilibrium 
polarization of the system in the absence of external excitation sources (/ = 0, g = 0). The 
angled brackets in Eqs. 10, J2J, and (0 denote the trace over reduced density matrix p(t), 
which is obtained by tracing out all environment degrees of freedom: (crx) = Tr(axp(t)), 
etc. In order to simplify the problem we assume the relaxation, Yz and dephasing, T, rates in 
(0-10 are time- independent, i.e. the rates are slowly varying functions on the scale of Rabi 
period which is of the order of l/f. 

Assuming that a high frequency driving amplitude F is sufficiently small we write the 
desired solution of Eqs. ©, and JSJ as: 

W = Z{t), (6) 



(oy) = Y(t) + Alt) cos(w t) + B(t) sin(w i), (7) 

{cr x ) = X(t) + C(t) cos(w t) + £>(t) sin(w t), (8) 

where X, Y, Z, A, B, C, and D are slowly varying, as compared to the high frequency u)q, 
quantities. 

As is known, a two-level system resonantly irradiated with a high frequency undergoes 
a low frequency Rabi oscillations. However, if the external low frequency excitation, G{t) 
is absent (g = 0), the Rabi oscillations are damped out. For this case we obtain at the 
degeneracy point (e = 0) the following solution: X — 0, Y = 0, C = B, D = —A, 

Z{t) = z* + z e~ Flt + e~ T2t (z x cosuj R t + z 2 sintu R t) , (9) 



A{t)=a*+a e r±t + e T ' 2t {a x cosu> R t + a 2 sinu> R t) , (10) 

B(t) = b* + b e- rit + er T2t (6i casu R t + b 2 sinu R t) , (11) 
where the relaxation rates, T\, T 2 , and the Rabi frequency uj r are determined by the equation 

(A - iT) 2 (A - iT z ) - (A - »r) f 2 - (A - iY z ) S 2 = 0. (12) 

Here A is the eigenvalue for the oscillation mode, e. g., Z{t) » e lXt . This cubic equation has 
simple analytic solutions only for two cases. Firstly, in the absence of a relaxation in Eqs. © 
- © Q?z = T = 0) we obtain Ti = T 2 = 0, uj r = \/S 2 + f 2 , where 5 is the high frequency 
detuning parameter, S = ujq — A e /h. Secondly, for 6 = we obtain Ti = T, T 2 = {T + Tz)/2, 
uj r = y/ f 2 — (r — r^) 2 /4. For the subsequent derivation we need only the quantities z*, a*, 
b* which can be determined as the steady state solution of the equations ©, (QJ, and JSJ): 

z* = P Z Z(0), a* - z*fT/(T 2 + S 2 ), b* = z*fS/(T 2 + S 2 ), where P z = rz yif 5 ^^ r - 

The quantity PzZq is the noncquilibrium polarization, i. e., the steady state difference 
of occupation probabilities between two energy levels in the case when the high frequency 
excitation is applied to the TLS. Therefore, in the absence of the low frequency excitation the 
Rabi oscillations decay with a rate T2 given by Eq. l|12f) . The main goal of our investigation 
is to obtain the persistent oscillations of the quantities Z(t) ©, Y(t) J7J, X(t) JSJ which can 
be detected with low frequency (compared to the gap) electronic circuitry. In what follows 
we show that a low frequency signal applied to TLS can sustain the persistent low frequency 
oscillations of the above mentioned quantities. The amplitude of these oscillations has a clear 
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resonance at the Rabi frequency with the width of the resonance being dependent on the 
damping rates of the system. 

The Eqs. ©,(01), and (J5J are analized in the presence of a low frequency excitation G(t). 
We insert Eqs. ©, 07 © into ©, 0, and and in accordance with the ideology of 
rotating wave approximation retain only the low frequency terms. We consider the force G 
to be small, and therefore neglect the terms which are of second and higher order in G. In 
addition, we keep only the terms which oscillate within the bandwidth of the Rabi frequency 
and neglect the terms which are of the order of Hf/A £ , HT/A £ , HT Z /A £ . Consequently we 
obtain the following set of equations for low frequency quantities: 

Z = fA- T Z Z (13) 



Y=^X-TY + 2gz*-^fB 

A F e 

X = -Y - TX + —fA 

h A 



(14) 
(15) 



A + 2TA + (0| + T 2 )A = 
f(T z -T)Z- T^gb* - S^ga* - ^gb* 



(16) 



B + 2TB + (5 2 + T 2 ) B = 
-SfZ + T^ga*-S^gb* + ^ga* (17) 

From these equations we can obtain the Fourier components for the slowly varying quantities 
A{u>) and B{uS), and for the low frequency persistent response of the TLS to a small low 
frequency excitation, Z(ui), Y(ui), and X(uu): C(uj) = B(ui), D(lu) — —A(u>), 

*»>-*">i^*^^5r (18) 

A( U ) = 9(^fP Z Z 0WTf - 2 [ - ^ *± (19) 

B(u) = -g{u)-f Pz Z j2 - r2 {p + {T + iuj)2) | j, (20) 

he-,. ( h \ 2 / . , „ „ s 



rH = a7K /am + (at) (*)^^4/ B H)( r + H, (21) 

X{w) = (£\ ±fA(u)(T + iw)- JL(2g(u)P z Z -lfB(uj), (22) 



where s(u) = (Sl B - to + iT)(fl R + u - iT)(iuj + T z ) - f 2 (T z — T); il R = ^S 2 + f 2 is the 
Rabi frequency in the absence of the damping (T z =7 = 0). ( : ) 

( 1 )In above expressions we disregard the Fourier components of the damping terms (Eqs. 1101 . 1111 ) 

since for the persistent signal g(t) they are unimportant. 
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Fig. 1 - Time evolution of < a z >. (thick) G=0, (thin) G/h = 1 MHz. The insert shows the 
undamped oscillations of < (j z > at G/h = 1 MHz. 



As can be concluded from these expressions the persistent low frequency oscillations of 
the spin components Z(t), X(t), Y(t) appear as the response to the low frequency external 
force, g(t), only in the presence of the high frequency excitation (/ ^ 0). Exactly at resonance 
(S = 0) A(u) — 0, and Z(ui) = 0, but B(ui) ^ 0. At this point the population of the two levels 
are equalized, and spin circularly rotates in the XY plane with frequency loq, with the center 
of the circle being precessed with the frequency lo of the low frequency external source. 

As an example we show below the time evolution of the quantity (az)(t) — Z(t), obtained 
from the numerical solution of the equations © , (@J , and 10 where we take a low frequency 
excitations as G(t) = Gcos^LOLt). The calculations have been performed with initial conditions 
(<r z )(0) = 1, (ct x )(0) = (ov)(0) = for the following set of the parameters: F/h = 36 MHz, 
A/h = 1 GHz, T/27T = 4 MHz, T z /2tt = 1 MHz, e/A = 1, Z = -1, <5/2tt = 6.366 MHz, 
lu l /Q r — 1. As is seen from Fig^ in the absence of low frequency signal (G — 0) the 
oscillations are damped out, while if G ^ the oscillations persist. 

The Fourier spectra of these signals are shown on FigEJ for different amplitudes of low 
frequency excitation. For G — the Rabi frequency is positioned at approximately 26.2 MHz, 
which is close to £Ir — 26.24 MHz. With the increase of G the peak becomes higher. It is 
worth noting the appearance of the peak at the second harmonic of Rabi frequency. This peak 
is due to the contribution of the terms on the order of G 2 which we omitted in our theoretical 
analysis. The comparison of analytical and numerical resonance curves calculated for low 
frequency amplitude, G/h = 1 MHz and different dephasing rates, T are shown on Fig|3] 
The curves at the figure are the peak-to-peak amplitudes of oscillations of Z(t) calculated 
from Eq. H18|) with g(u>) — g{5(ui + w^) + 5(ui — ui^))/2, where 5(ui) is Dirac delta function. 
The point symbols are found from numerical solution of Eqs. (0 , (@J , (JSJ . The widths of the 
curves depend on T (see the insert) and the positions of the resonances coincide with the Rabi 
frequency. A good agreement between numerics and Eq. ^3 as shown at Fig. |3 is observed 
only for relative small low frequency amplitude G/h, for which our linear response theory is 
valid. 

The key point of the method we described above is that it allows for the detection of the 
high frequency response of a TLS at a frequency which is much less than the gap frequency. 
The low frequency dynamics of the quantities (ax) , (fr), (<Jz) bears the information about 
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Fig. 2 - Fast Fourier transform of < o z > at different amplitudes G/h of low-frequency field. 

the relaxation, Tz, and dephasing, T, rates. 

The experimental realization of this method (the detection scheme) depends on the prob- 
lem under investigation. For example, the method can easily be adapted for NMR. In this case 
the quantities (cr z ) , (ffx), (°y) arc the longitudinal, Mz, and the transversal polarizations 
of the sample, Mx, and My. Indeed, from the comparison of the Hamiltonian @ and the 
equations of motion |(3}, Q, (J5J with those for nuclear spin (H = —ftB, dfl/dt = 7 [/IB], where 
7 is the gyromagnetic ratio), it is clear that the NMR case corresponds to the polarization 
of a sample with a field B along the z axes with a high frequency excitation Bicos(ujQt) 
and a low frequency probe G(t) being applied in the ZX plane at an angle 8 to the z axis: 
B z = B + coseiBxeosiuot) + G(t)), B x = -sind(B lC os{u t) + G(t)), B Y = 0. This anal- 
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Fig. 3 - Peak-to-peak amplitude of < a z > pp undamped harmonic signal at t « 700 ns vs ul/Qr- 
G/h = 1 MHz, r z /27r = 1 MHz, (solid curve) r/27r = 6 MHz, (dashed curve) r/2vr = 4 MHz, (dotted 
curve) r/2-7r = 3 MHz. The curves and the width of resonance, AW at 1/V2 level (see the insert) 
are calculated from Eq. 1181 : the symbols are found from numerical solution of Eqs. 
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ogy allows for the direct application of Eqs. ^18 fl -l]22 p to the nuclear spin with the following 
substitutions: A e = jhBo, e/A e = cosO, A/A e = sin9, f = jBisinO, g(uj) = jG(ui)sin6, 
e/A = cot 9, Zq = Mo - the equilibrium magnetization of a sample. In NMR all low frequency 
components of the magnetization, Z(oj) 1|18|) . Y(lu) l|21l) . X(lo) (|22() . and their combinations 
are accessible for the measurements. 

Our method can be directly applied to a persistent current qubit, which is a superconduct- 
ing loop interrupted by three Josephson junctions [16], [17]. For these qubits the successful 
experimental implementation of low frequency readout electronics has been demonstrated 
[18-20]. The average current in the qubit loop is proportional to the low frequency part of the 
quantity (az)(t), which is directly connected to the probabilities of occupation of the ground, 
P-(t), and the excited, P+(t) states: (az)(t)—P-(t)-P + (t). Therefore, this current can be 
detected through the variation of its magnetic flux either by a DC SQUID [21] or by a high 
quality resonant tank circuit inductively coupled to the qubit [18-20]. 

In conclusion, we proposed a method to study the Rabi oscillations in a dissipative TLS 
by irradiating it simultaneously with high, resonant, and low frequency. The low frequency 
response of a system has a clear resonance at the Rabi frequency with the resonance width 
being dependent on the damping rates. Therefore, the method allows for the experimental de- 
termination, in the low frequency domain, of the relaxation and dephasing rates of dissipative 
two-level systems. 
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